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Abstract- The random medium model is used to interpret the 
polarimetric active and passive measurements of saline t ice. The 
£lLr is described as a host ice medium embedded with 
randomly distributed inhomogeneities and the ^Jrl^ng sea 
•.ter is considered as a homogeneous half-space. The scattere 
ta Uie ice layer are modeled with an ellipsoidal correlation 
function The orientation of the scalterers is vertically aligned 
S azimuthally random. The strong permittivity fluctuation 
theory is employed to calculate the effective permittivity and the 
distorted Born approximation is used to obtain the polanmetnc 
scattering coefficients. We also calculate the thermal emissions 
Sd on the reciprocity and energy conservation pnnap'esjhe 
effects of the random roughness at the air ice, an ice- 
interfaces are accounted for by adding the surface scattering to 
. v volume scattering return incoherently. The above theoretical 
model which has been successfully applied to analyze the radar 
backscatter data of the first-year sea ice near Point Barrow 
AK, is used to interpret’ the measurements performed in the 
CRRELEX program. 

INTRODUCTION 

Understanding electromagnetic wave interaction with saline ice 
provides physical insight for interpretation of active and passive 
remote sensing data of sea ice. Various theoretical models have 
been developed to characterize the scattenng and the emission 
properties of sea ice [1]. In this paper, the random medium 
model [2] is used to interpret the measurement data collected 
at the US Army Cold Regions Research and Engineering Lab- 
oratory (CRREL) under the CRRELEX program. To simulate 
sea ice, saline ice was grown under partially controlled condi- 
tion at CRREL and active and passive microwave signatures 
of the saline ice were measured together with ice characteristic 
parameters. 

Under Quiescent condition, the background ice grows in 
columnar fonn and saline water is trapped between ice platelets 
in the form of brine inclusions which are usually ellipsoidal. The 
ice tends to grow vertically downward rendering the ellipsoidal 
inclusions aligned preferably in the vertical direction and the 
crystallographic c-axes parallel to the Horizontal plane, n t .s 
case, the c-axes are. however, random m the horizontal direction. 
The saline ice laser is modeled as a random medium containing 
ellipsoidal scattered vertically aligned but randomly oriented 
in the horizontal direction. The strong fluctuation t ieor\ is 
employed to calculate the effective permittivity and the distorted 
Bom approximation is used to obtain the polarimetric scattenng 
coefficients. The emissivity of a layer of sea ice overlying sea 
water can be obtained through the principles of reciprocity 
end energy conservation [3] . In other words, the emissivity 
is obtained as one minus the reflectivity. To account for the 

L "’in method 

1 ice-water 


interfaces, and add incoherently to the volume scattering return 
obtained with the distorted Bom approximation. This model 
has been successfully applied to study the radar backscatter for 
the thick firet-year sea ice data measured near Point Barrow, 
Alaska [5j. In this study, we further employ this model to 
interpret the multi-frequency active and passive measurements 
taken during the CRRELEX experiment. 
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Figure 1 Scattering configuration. 


model description 

Multi-Layered Random Medium Model 

The fully polarimetric multi-layered random medium model 
for sea ice has been developed with the wave theory. The 
layer configuration is shown in Fig. 1. Region 0 represents 
the homogeneous air above sea ice. The sea ice m region 
contains vertically aligned ellipsoidal scalterers such as sa-ine 
ice with brine pockets or air bubbles which has spatially varying 
permittivity <i(f). Region 2 is the underlying sea water w.th 
homogeneous permittivity €2 All three regions are assume o 
have the frec-space permeability fi 0 . 

An exponential correlation function locally corresponding to 
a scatterer is used in the derivation of the effective permittivity 
with the strong fluctuation theory extended to account for 
the ellipsoidal shape and the horizontally random orientation 
of the scalterers. The local correlation function of the brine 
inclusion has three different correlation lengths related to tie 
three different axial lengths of the ellipsoid as shown m Mg 
2 The azimuthal orientation of the scatterer is depicted locally 
with the local z'-axi; parallel to the global z-axis and an Eulenan 
fluctuation angle 6t corresponding to the angle between the c- 
axis and the global horizontal axis (x-axis). When the average 
process is performed over the random orientation angle Of, 
the effective permittivity is a uniaxial tensor with vertical 
optic axis describing the effective azimuthal symmetry ol the 
inhomogeneous sea ice. When a single set of correlation lengths 
is used, for simplicity, to model the brine inclusions with vanous 
sizes, these correlation lengths are considered as related to the 
effective size of the scalterers. After the effective permittivity 
of the sea ice is calculated, the complete set of polanmetnc 
scattering coefficients are obtained under the distorted Bom 
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approximation with the anisotropic dyadic Green’s function for 
the layer media. From the scattering coefficients, the Mueller 
matrix is computed to plot the polarization signatures of sea 
icc. 

Under the distorted Bom approximation the polarimetric 
scattering coefficients, obtained as the correlation of the scat* 
tered field, takes the form 

(E 0 s (r)T 0s <?)) = 

^ k* [ dor p(<?f) I dfi[ dr°C^ jk , m (r u r° v .oj) 

. “ JO Jv l Jv \ 

( 1 ) 


l 



where the uniaxial effective permittivity i e jj\ for region 1 is 
used to calculate the mean dyadic Green’s functions and the 
mean fields (F](f)). The correlation function for region 1 
in (!) is the jklm-th element of fourth-rank 

correlation tensor Cci(n ,r^; Qj) defined as 

Q ljWm (fi,7t;?/) = (fljk(n) |^/(n)) (2) 

In (]), the probability density function of scatterer orientation 
pW>/) is simply 

?(<*>/) “ TjZ ^ 

for complete random orientation horizon tally with no preference 
in azimuthal direction, and k 0 = ujyjvoto Is the free-space 
wavenumber with angular frequency ur. 


Effective Permittivities 


The strong fluctuation theory is used to derive the effective 
permittivities of the ice layer. The singularities of^the dyadic 
Green’s functions in the bilocal approximated Dyson s equations 
are" accounted for; and the low-frequency approximation is 
applied to obtain the results for the anisotropic sea ice media. 
The derivations are done in the frequency domain with the 
Fourier-transform method. 

The ice layer is described with a normalized local correlation 
function of the form 


i(rO = exp - 


x' 2 
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with correlation length £ x / t £^f, and in the local coordinates 
corresponding to the minor, the meridian, and the major axes of 
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the scatterer. In this model, the correlation lengths are related 
to the effective size and shape of the scatterers. The effective 
permittivity of the ice layer is composed of a quasi-static part 
and a scattering part which accounts for the modification in the 
wave speed and attenuation 

Ce//1 ” sc//l (^l)] (5) 

The effective dyadic scatterer \ t jj\ is obtained from the volume 
integration of the dyadic Green’s function and the correlation 
function, the auxiliary permittivity 7 g \ and dyadic coefficient 
S § 1 characterizing the source singularity of the dyadic Green’s 
function are determined by the condition or secular term 
elimination. 


Scattering Coefficients 

For polarimetric backscattering, the scattering coefficients 
are defined by 
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where A is the illuminated area, subscript i is for the incident 
wave, s represents the scattered wave, and /i, v , r, and * 
can be h or v which respectively stands for the horizontal or 
vertical components of the electric fields. The components of 
the scattered field in (6) are obtained by measuring the h and 
the v returns while the incident field is transmitted exclusively 
with h or v polarization. For backscattering from reciprocal 
media, the h return obtained by transmitting v is the same as 
the u return obtained by transmitting h. 


Rough Surface Contribution 

To account for the effects of rough interfaces at z = 0 
and z = -d\, we first approximate the surface statistics by 
a Gaussian correlation function with standard deviation o and 
correlation length £, then employ the small perturbation method 
to obtain the first-order backscattering field From the interfaces. 
For example, the surface scattering from the air-ice interface is 
given by 


°hh = cos 4 8, 

a vv = 4k^cT'(~ cos 4 8, 



-klc-sxr? 8, 


(*? - *?)(*? sin 2 Oj + kj.) 
(*?*, + * 5*u) 2 



(7a) 

(76) 


where 0, and Jq denote the angle of incidence and the wavenum- 
ber for the ice layer, respectively. The surface scattering is inco- 
herently added to the volume scattering return calculated based 
on the method described previously. 


Emissivities 

The Kirchhoff law is invoked to obtain the emissivity of 
ice under isothermal condition from the active bistatic scattering 
coefficients. 


e/i(0. 4) = 


ev{0,<j>) = 



(50) 
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where Thh ^’ represents the bistatic scattering coeffi- 
cient for an h-polarized incident wave propagating in the direc- 
tion described by (0,<$) and scattered into the {S',#) direction 
with 5-polarization. 

DATA INTERPRETATION 

The theoretical results from the random medium model are com- 
pared to the CRRELEX data. For the case under consideration, 
region 0 is air, region 1 is the sea ice layer containing ellipsoidal 
brine inclusions, and region 2 is sea water. The thickness of 
the ice layer is obtained by ground truth measurement. The 
fractional volume of the brine inclusions is calculated from the 
bulk ice salinity and temperature. Results from publications are 
used to estimate the frequency-dependent permittivities of the 
ice background, the brine inclusions, and the underlying sea wa- 
ter. The theoretical calculations performed with and without 
the rough interfaces contributions are compared. 

The local correlation lengths, which are not defined in the 
same manner as the global correlation lengths in the Perovich's 
estimation, corresponding to the effective size of the scatterers 
need to be determined. In our preliminary calculations, instead 
of adjusting the correlation lengths to match the experimental 
data, they are chosen from the characterization of the first-year 
sea ice obtained by a comparison with the 9 GHz measurements 
made at Point Barrow, Alaska [2]. This approach, based on the 
theoretical model, also serves the purpose of referring the remote 
sensing signatures of the saline ice artificially grown at CRREL 
to the sea ice in the natural conditions. 

The theoretically calculated emissivities are in good agree- 
ment with the data. This indicates that the estimated effective 
permittivity tensor, which affects the emissivity even in the ze- 
roth order, are close to that for natural sea ice. At lower fre- 
quencies where the electromagnetic wave can penetrate the ice, 
the emissivities can be sensitive to the ice thickness. In this 
case, the h polarization is more sensitive to the thickness than 
the v polarization which usually has higher attenuation. For 
information pertaining more to the first or higher order effects, 
active signatures should be considered. 

The backscauering coefficients generated theoretically have 
also been compared with that observed in the CRR.ELEX ex- 
periment for co-polarized and cross-polarized returns only, since 
measurements of other elements in the covariance matrix are not 
available. At 5 GHz. the theoretical results overestimate the CR- 
RELEX data at larger incident angles. This indicates that the 
local correlation lengths of the grown sea ice probably smaller 
than those observed at Point Barrow. At smaller incident an- 
gles. the measurements are higher than the theoretical values 
without rough surface contributions. At 10 GHz. the results, 
calculated using the above correlation lengths, are all higher 
tr.an the measured values provided to us for the backscauering 
* r om the CRRELEX experiment. 
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Scattering of Electromagnetic Waves from a Dense Medium 
Consisting of Correlated Mie Scatterers with Size Distributions 
and Applications to Dry Snow 

L. Tsang 

Electromagnetics and Remote Sensing Laboratory 
Department of Electrical Engineering 
University of Washington 
Seattle, WA 98195, USA 

J. A. Kong 

Department of Electrical Engineering and Computer Science 
Massachusetts Institute of Technology 
Cambridge, MA 02139, USA 

Abstract- The scattering of a plane electromagnetic wave obliquely incident on a layer 
of dense medium consisting of dielectric spherical particles of finite sizes and with size 
distributions is studied- The spherical particles are of sizes comparable to wavelength 
so that Mie scattering is used to describe the single particle scattering characteristics. 
The coherent wave is studied with quasi crystalline approximation using the cross pair 
distribution functions of multiple sizes governed by Percus-Yevick approximation. The 
incoherent scattered wave is calculated with the distorted Born approximation with the 
result expressed in terms of a product of the T-matrices of particles of different sizes and 
permittivities and the Fourier transform of the cross pair distribution functions. The co- 
herent wave effective propagation constants, the attenuation rates and the backscattering 
coefficients are illustrated numerically, with examples chosen to illustrate microwave and 
millimeter wave scattering from snow cover in the frequency range of 5 GHz to 95 GHz, 
and mean grain radius between 0.03 cm to 0.09 cm. Salient features of the numerical 
results for scattering from snow with size distribution are: 1. Correlated dense medium 
scattering is less than independent scattering at low frequency, a fact that is consistent 
with controlled laboratory experiment. 2. Scattering from dense medium of a Rayleigh 
size distribution with an average radius can be much larger than the case of monodisperse 
particles of the sizes identical to that average radius. 3. The scattering attenuation rate 
increases rapidly with frequency at low frequency regime and begins to level off at high 
frequency regime. 4. The coherent wave scattering attenuation rate is large in snow at fre- 
quencies above 15 GHz Indicating large optical thickness and the importance of multiple 
scattering. Comparisons with extinction measurements of dry snow at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz are made. 


L INTRODUCTION 

Propagation and scattering in dense media has applications in remote sensing of 
geophysical terrain and nondestructive evaluation of composite materials [1-13]. 
In a dense medium, particles occupy an appreciable fractional volume. A main fea- 
ture of a medium with densely packed scatterers is that Foldy’s approximation [14] 
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which has been successful in propagation problems through sparse media, is not 
applicable. It has been shown theoretically [1-2,4-6,9,15-16] and through labo- 
ratory experiments [3] that the effects of correlated scattering between scatterers 
must be taken into account. 

Propagation and scattering in dense media has been studied with the quasicrys- 
talline approximation [1-2, 4-6], and the quasicrystalline approximation with co- 
herent potential [6-8] for the first moment of the field and the correlated ladder 
approximation for the second moment of the field [7,9,13]. The dense medium 
theory is developed for media with particles of finite sizes. Because of the finite 
sizes of the particles, their positions are correlated. The dense medium radiative 
transfer theory has also been developed from these approximations to study mul- 
tiple scattering effects in the incoherent wave. Some of the studies of correlated 
scattering have been limited to the case of low frequency [7,9] and/or particles 
of identical sizes [1-2,4-6,10-11]. However, to approximate physical reality in 
geophysical terrain scattering, size distributions of particles must be taken into 
account. The case of monodisperse particles and polydisperse scattering generally 
can give different results for independent scattering [14,17] as well as correlated 
scattering [8,16]. High frequency results are also important for example at fre- 
quencies above 15 GHz for snow cover. 

In this paper, we study the scattering of a plane electromagnetic wave obliquely 
incident on a layer of dense medium consisting of dielectric spherical particles with 
size distributions. The spherical particles are of sizes comparable to wavelength so 
that Mie scattering is used to describe the single particle scattering characteristics. 
Thus the theory is applicable to particles of arbitrary size. The coherent wave is 
studied with the quasicrystalline approximation using the cross pair distribution 
functions of multiple sizes governed by Percus-Yevick approximation [18-20]. The 
integral equation of the quasicrystalline approximation gives rise to the generalized 
Lorentz-Lorentz law and the generalized Ewald-Oseen extinction theorem [4,21]. 
The generalized Lorentz-Lorentz law consists of a homogeneous system of equa- 
tions. The resulting dispersion relation is an equation governing the effective prop- 
agation constant for the coherent wave. The generalized Ewald-Oseen extinction 
theorem is an inhomogeneous equation that relates the transmitted coherent field 
to the amplitude, polarization, and direction of propagation of the incident field. 
The name “generalized” denotes the fact that the classical Lorentz-Lorentz law 
and the Ewald-Oseen extinction theorem were derived for low frequency effective 
permittivity for the total field with the imaginary part accounting for absorption 
loss only. In the present “generalized” case, the effective propagation constant is 
derived for high frequency when the particle sizes are comparable to wavelength. 
Also the imaginary part of the effective propagation constant includes attenuation 
due to scattering diverted to other directions as well as absorption. 

The calculations of the cross pair distribution functions are based on a spe- 
cial case of Percus-Yevick approximation that has no inter-particle forces except 
for non-interpenetration [8,17-19,22]. The Percus-Yevick cross pair distribution 
functions of multiple sizes have recently been shown to compare very well with 
Monte-Carlo simulations [23]. The incoherent scattered wave is calculated with 
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the distorted Bom approximation with the result expressed in terms of a product 
of the T-matrices of particles of different sizes and permittivities and the Fourier 
transform of the cross pair distribution functions. Numerical examples are chosen 
to illustrate microwave and millimeter wave scattering from snow cover in the 
frequency range of 5 GHz to 95 GHz and mean grain radius ranging from 0.03 cm 
to 0.09 cm [24]. The Rayleigh size distribution and the modified gamma size dis- 
tribution are used to illustrate the scattering. Salient features of the numerical 
results for scattering from snow with Rayleigh size distribution are: 1. Correlated 
dense medium scattering is less than independent scattering at low frequency 
that is consistent with controlled laboratory experiment [3]. 2. Scattering from 
dense medium of a Rayleigh size distribution with an average radius can be much 
larger than the case of monodisperse particles of the sizes identical to that av- 
erage radius. Particles with larger size can contribute to scattering even though 
they are fewer in number. 3. The scattering attenuation rate increases rapidly 
with frequency at low frequency regime and begins to level off at high frequency 
regime. This is to be contrasted with Rayleigh scattering where scattering at- 
tenuation just increases rapidly with frequency. 4. The coherent wave scattering 
attenuation rate can be large at frequencies above 15 GHz based on the Rayleigh 
size distribution indicating large optical thickness and the importance of multiple 
scattering in snow at high frequencies. Comparisons are made with extinction 
measurements of dry snow at 18 GHz, 35 GHz, 60 GHz, and 90 GHz. Good agree- 
ment is obtained and the dense media theory is shown to be able to account for 
the frequency dependence of extinction of dry snow from 18 GHz to 90 GHz. 


II. FORMULATION 

Consider a plane electromagnetic plane wave impinging in direction (n - 0,, 4>i) 
upon a half space of spherical discrete scatterers of multiple sizes and multiple 
permittivities (Fig. 1). Let sj be the species index j = 1, 2, . . . , L. Each species is 
described by the permittivity c s j, radius a s j and number density n s j. The num- 
ber density n s j is related to the volume fraction of species f s j, f s j = 

The scatterers are embedded in the lower half space (region 1) with background 
permittivity equal to e which is the same permittivity as that of the region 0 

(Fig. 1). 


The electric field of the incident wave is given by 

Ei{r ) = (Ej id + E hi 4>i) (1) 

where subscripts v and h denote vertical and horizontal polarizations respec- 
tively, denotes the downward propagating incident wave vector and and 
are the associated polarization vectors 

kid — &(sin #i cos <f>% £ + sin s i n <t>i V — cos (2) 

= % d* kiy y ki z 2 

$ id = -i cos cos <j>i - y cos 6{ sin fc - z sin 0, (3) 

sin <f>j + y cos <f>{ (4) 
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In the following, our notation follows that of [6], pp. 446-454 and pp. 506-517. 
Under the quasicrystalline approximation for multiple species, the coherent field 
obeys the equation 


t D^ s, \ri) = n *j / 9»M r j ~ 7 l) °( k7 l r j) 

Mj = 1 JV » 

+ &*** 



( 5 ) 


where uK*')(r/) is the conditional average of the exciting field of particles at r/ of 

species s/ conditioning on the particle at r/, T * is the matrix representation of 
the T matrix of species sj and F(Jfcr/rj) is a transformation matrix transforming 
from vector spherical waves with center at Tj to vector spherical waves with cen- 
ter at r/ as given on p. 452 of [6]. It is an exact wave transformation applicable 
to near field, intermediate field, and far field regions. The function gspfij - n) 
is the cross pair distribution function of two species Sj and s/. These can be 
computed readily from the Percus-Yevick approximation. These were originally 
derived for molecular dynamics of fluid mixtures. The result that we use is the 
special case when intermolecular forces are zero, and molecules cannot interpen- 
etrate each other. In (5) the volume of integration Vji for drj is the lower half 
space excluding a “hole” of radius ( a s ■ + a*,) around vector r/. This is because 

9tjs,(rj - r/) = 0, for | \rj - r/| < a S} + a s , (6) 



Figure 1. Geometry of the problem: an incident electromagnetic plane 
wave impinging upon a dense medium consisting of spherical 
particles with a size distribution. 
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In (5) a,„ c is the column matrix representing the coefficients of the incident 
wave when expanded into spherical waves, as given on pp. 507 and 514 of [6] . 

To solve (5), we assume the solution 

tzK")(n) = (7) 

where is a column matrix representing coefficients of the exciting field and 

~K 4 is the downward going effective wave vector 

K d = K x x + K y y-K z z (8) 

and the effective wavenumber is 

K = Jk* + Kf + Kl (9) 

Substitution of (7) into the integrand of the integral (5), integrations can be 
performed [4,6]. The integral yields two terms that correspond to two waves with 
dependencies 

(io) 

811(1 c » K x xi+iK vVl -i(k 2 -Kl-Kl) l f 2 z, ( 11 ) 

The first wave given by (10) propagates with K d and is to be balanced with the 
term w s, ( 7 l) on the left hand side of elation since the y both have the same 
propagation constant. Balancing these two terms give rise to the generalized 
Lorentz-Lorentz law which can be used to calculate the effective propagation con- 
stant K. The second wave of (11) with incident wave wavenumber k in the square 
root will cancel the second term on the right hand side of (5) which is the incident 
wave term. This cancellation or extinction is the generalized Ewald-Oseen extinc- 
tion theorem [6,25]. For the two terms to cancel out, the wave vector of (11) must 
exactly balance the incident wave vector. This gives rise to the phase matching 


condition 



■ie* 

II 


* 

II 

so that from (8), 

(9), (10), we obtain Snell' 

constant, 



K smdt = ksin8i 0^) 

K z = (/v 2 — fc 2 sin 2 0,) 1/2 (15) 

We call both laws generalized because the concept has been generalized to coherent 
wave at high frequency while the classical laws were derived for low frequency case. 

A. Generalized Lorentz-Lorentz Law 

The column matrix for the exciting field can be expressed in terms of 

Cj 

y-(*j)(M) y n (*>)( JV ) coefficientsjvhere (M) stands for M type of vector 

spherical wave and ( N ) stands for N type of vector spherical wave and n is the 
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multipole index. The y n ( ^ )(M) and V„ ( * i)(A ° coefficients satisfy a homogeneous 
set of equations. By setting the determinant equal to zero, one gets a nonlin- 
ear algebraic equation for the effective propagation K. Because the medium is 
statistical isotropic, we have the same homogeneous set of equation for arbitrary 
polarization and incident angle of the incident wave. The homogeneous set of 
equations is 

= —2* £ n„ 2n + l)S p (i, W;.,) 

*/=l »».P 

{a(l,n| - l,u\p)A(n,u,p)Ti t,){M) Y^ M)it,) 

+ a (1, n| - 1, u\p,p - 1) B(n, «/,p) 4 $l)(N) 

= _ 2 „ £ ^ £ (2ri + l)^ p (A:,/s:|^ i , l ) 

x {a (1, n| - 1, v\p,p - 1) B(n, u,p) r<*' )(M) r„ ( *' )(M) 

+ « (1, n| - 1, u\ p)A(n, u,p) T<" )(N) Y ^ ,)(N) } 

where o(l,n| - l,u\ p) and a(l,n| - l.Hp.P" 1) are in terms of Wigner 3-j 
symbols and can be found on pp. 449-450 of [6]. In (16), (17) 

, i/(i / + 1) + n(n + 1) — p(p + 1) /.p, 

A (n, i/,p) K 19) 


(16) 


(17) 


v x/(n + j / + p + 1)(»/ +p - n)(n + p - »^)(n + v - p + 1) /1ft ^ 

B(n, v,p) = — 

Sp (it, K\R t jt t ) — ~ [^(fc-R*^) jp^'-R^*,) “ Khp(kR SjSl ) jpiKRsjS , )] 

+ I dr r 2 [ff<ji,(r) - l] h p (kr)j p (Kr) 
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For each species sj, we terminate at a multipole N Ct depending on the radius 
a of that species. Thus (12) and (13) form N e homogeneous equations for 
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the N e unknowns y^*'^ \ n = 1 , 2, . . . , N $l and s/ = 1,2, . . . ,L. 

Thus N e = EiLi 2JV*,- . . 

The generalized Lorentz-Lorentz law determines the effective propagation K. 
Setting the determinant of (16)-(17) equal to zero gives the nonlinear equation for 
the effective propagation constant K. After K is calculated, then the complex 
transmitted angle 9i and Kz can be determined from (12)— (15). 

The N e coefficients y n ( * ,)(M) and Y^‘ m are then reduced to only one arbi- 
trary constant that is to be determined by the incident wave. That one equation 
to determine the arbitrary constant is provided by the generalized Ewald-Oseen 
extinction theorem which gives a single inhomogeneous equation. 

B. Generalized Ewald-Oseen Extinction Theorem 

The generalized Ewald-Oseen extinction theorem is obtained by balancing the 
incident wave term of the second term of (5) and the term of (11) that is a result 
of the integral in (5). The result matrix equation can be reduced to one single 

equation for the set of coefficients and Y^ S ‘ ){N] . The equation depends 

on whether the incident wave is vertically polarized or horizontally polarized. 


1. Vertical polarized incidence 

In this case E vi / 0 and E hi = 0. The coherent transmitted wave is also 
vertically polarized. The generalized Ewald-Oseen extinction theorem gives rise 
to the following inhomogeneous equation 


. (2n + 1) f T ( < ,KM) r ( < ,)(M)(V) P n 1 (cos(g, - fit)) 

,1 ® ,, n(n + 1) L n " |sin(0j — $t)| 

8 (— 1 n 

f s(e , _ + n( „ + !)/>„(»<*, - »,»}] 


(Kz ~ ~ kjz ) kj z k 

27 r 


(24) 


where the additional superscript (V) denotes the vertically polarized case. Note 

that the set of coefficients yM^^ and y^ N ^ V \ n = s, = 

1 L have been determined to within one arbitrary constant by the generalized 
Lorentz-Lorentz law of (16) and (17). Thus (24) provide the last equation that 
determines all those coefficients uniquely. 


2. Horizontal polarized incidence 

For this case E vi = 0 and E hi 4 0. The coherent transmitted wave is hori- 
zontally polarized. The inhomogeneous equation for the generalized Ewald-Oseen 
extinction theorem is 


EE"- 


,(2n±l) 


n(n + 1) 


r (*i)(M) v (5,)(M)(ff) 

‘n J n 
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where the additional superscript (/f) denotes horizontal polarization. Thus the 
homogeneous system of equations of (16), (17) combined with the inhomogeneous 

equation (25) determines the coefficients and y^r'K#) uniquely 

for the use of horizontally polarized incidence. 

The case of an incident wave of arbitrary polarization can be calculated by 
linearly combining the cases of vertical and horizontal polarized incidence. 


C. Bistatic Scattering Coefficients for Incoherent Field 

The incoherent field can be studied with the distorted Born approximation 
or the first order smoothing approximation [4,6]. From the incoherent field, the 
bistatic- scattering coefficients for the incoherent field can be calculated. The 
bistatic scattering coefficients are defined by [6] 


r \ 47rr 2 |Z r ^| 2 

7 = |f£|2 AoCOS< ?. 


(26) 


where £*p is the incoherent scattered field of polarization /?, E& is the incident 
field of polarization a, A§ is the area of target and r is the distance of observation 
point from target. In the following, we shall give the result for bistatic scattering 
in the plane of incidence <f>s = 4>i + tt and 0 < 6$ < tt/ 2. Because of statistical 
isotropy of random distribution of spherical scatterers, bistatic scattering in the 
plane of incidence do not have depolarization for vertically polarized incidence or 
horizontally polarized incidence. Under the distorted Born approximation, 7t>v 
and 7 m are given as follows; for 0 < B $ < tt/2 

7v*(0i, 4>S = K + 4>i\ 9i, = 2K" cos 8 \ W ‘PV‘M 

*i * * 


+ j^2 cos Q 
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w {V) (a Vr n (2n + 1) f («,)(A/)(s,)(M)(V) P,j(cos(fl, - 6 t )) 
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and 

U S}Sl (9»,4>» = * + <t>i\8i,4>i) = 2 Im^i fl) S-oo ^ h jS,(F) " l ] e ' { ** 7 ‘ d ~ T ' ) f 


= 2Rfe)^dP=-^ + I,) 


( 31 ) 


where superscript double prime denotes imaginary part. In (31), Tfy is given 
by (8), ( 12)— (15) and 

k s ~ Jt(sin 6 S cos <f>$x + sin 0 S sin <f> s y + cos 8$z) (32) 

is the wave vector of the scattered wave in region 0. In (31) Sj , s/ = 

1,2, is the Fourier transform of the cross pair distribution function g SjSJ 

. . . t m i . . r 


minus 1. It is also called structure factor. 

too 


= (2^)3 - 1 1 < 


->?r 


(33) 


Thus the bistatic scattering coefficient is proportional to the Fourier transform 
of the cross pair distribution function, a concept that has been used extensively 
in the investigation of molecular structures of liquids and amorphous solids by 
X-rays [22]. 

For multiple sizes of spherical particles, the cross pair distribution function gs^si 
are isotropic, so that gs^s^) depends on |r| only. Hence, the cross structure 
factor fl* j5/ (p) also depends on |p| only. It is useful to note that under the 
Percus-Yevick approximation, closed from expressions are available [18,19,22] for 
Jf*j*/(p)i s l “ 1) 2, * . * , L. 

In (27) and (28) 


F„=(l 


(34) 


where K” = In^if*). 

We have inserted in (34) the thickness d of the layer containing random discrete 
scatterers. The results of (34) only include volume scattering and ignores any 
possible reflections from the bottom boundary. 
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III. NUMERICAL RESULTS AND DISCUSSION 


In this section, we shall illustrate the numerical results of the effective propagation 
constants and backseat tering coefficients using parameters of snow. The results 
are illustrated with the Rayleigh size distribution which is a special case of the 
modified gamma size distribution. The size distribution is represented by n(a), 
which is the number of particles per unit volume with radii between a and a + da. 


It is given by 



(35) 


where (a) is the mean 
defined by 


radius. The fractional volume 

f - J 0 da y ° 3n ( a ) 


and the mean radius are 


(36) 


J n °° daan(a) 

W ~ Jq° dan(a) 


(37) 



Thus, given the fractional volume /, the Rayleigh size distribution contains only 
one parameter, viz., the mean radius. In (38) a m is the mode radius where the 
size distribution is at a maximum [20], and as indicated in (38), is approximately 
equal to 0.798(a). In our numerical calculations of Figs. 2 to 11, we discretize 
the Rayleigh size distribution into 20 different equally spaced sizes (i.e. L = 20 
species). The pair distribution functions for a medium with 20 sizes are then 
calculated by the Percus-Yevick approximation [20]. The maximum size for these 
20 sizes is chosen such that at least 99.7% of the scattering contributions and 
the fractional volume of particles are included. That usually occurs around 2.5 
to 3 times the mean radius. Comparison will be made between the scattering 
of a medium with Rayleigh size distribution of a specified value of (a) and the 
case of monodisperse particles with radii equal to that same value. The results 
of QCA will also be compared with that of independent scattering. For the case 
of independent scattering, the extinction rate is 


(*e independent = 


£n 4> £(2n + l)Re 
tj =l 


(m)s 3 (n)»; 


(39) 


In Figs. 2 and 3, we plot respectively the effective loss tangent 2K(/K T and the 
normalized phase velocity k/K r as a function of fractional volume at frequency = 
37 GHz for two Rayleigh size distributions with (a) = 0.05 cm and (a) = 0.09 cm. 
We note that the effective loss tangents first increase with fractional volume, rise 
to maxima and then decrease with further increase in fractional volume. These 
are consistent with controlled optical experiments [3] and controlled microwave 
experiments [26]. The peak of the larger size of (a) = 0.09 cm occurs at / = 0.19 
and the peak for (a) = 0.05 cm occurs at / = 0.16. This indicates that scattering 
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by larger particles are “more independent.” “More independent” means that for 
scattering by large particles, the “hole correction” for distribution functions very 
often are sufficient to characterize the scattering [ 6 ]. Figure 3 shows that the phase 
velocity decreases with fractional volume with smaller velocity for larger particles 
indicating that scattering decreases phase velocity. In Fig. 4, the back scattering 
coefficients are plotted as a function of fractional volume for a layer thickness of 
d = 3 cm. The backscattering coefficients are = cosd, -ypp with 0 = v,h. A 
small thickness is chosen because the model of incoherent scattering is limited to 
distorted Born approximation which is valid only for small optical thickness. The 
associated extinction rates of Fig. 2 show that multiple scattering of incoherent 
waves become important for thickness larger than 3 cm as the optical thickness 
becomes large. The results of Fig. 4 indicate that backscattering coefficient can 
be as large as —4 dB for a small layer thickness of d — 3 cm. 



Figure 2 . Effective loss tangent 2 Ki/K r of QCA as a function of frac- 
tional volume, frequency is at 37 GHz. The permittivity of 
the particles is e* = (3.2 + i0.001)eo- Two cases of Rayleigh 
size distributions are considered with (a) = 0.05 cm and (a) = 
0.09 cm. 

■mi redlM l »a>. Ml a*. 



Figure 3 . Normalized phase velocity k/Kr of QCA as a function of frac- 
tional volume. Frequency is at 37 GHz. The permittivity of 
the particles is €$ = (3.2 + tO.OOl) 60 . Two cases of Rayleigh 
size distributions are considered with (a) = 0.05 cm and (a) = 
0.09 cm. 
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Figure 4. Backscattering coefficient in dB of QCA as a function of 
fractional volume, frequency is at 37 GHz. The permittivity 
of the particles is e t = (3.2 + t'0.001)e 0 . Two cases of Rayleigh 
size distributions are considered with (a) = 0.05 cm and (a) = 
0.09 cm. The layer thickness is d = 3 cm and incident angle is 

30°. 

In Figs. 5, 6, and 7, we show the coherent wave attenuation rates for a smaller 
mean grain radius of 0.035 cm. The corresponding mode radius is 0.0279 cm. 
Figure 5 compares the extinction rate 2A'j between independent scattering and 
QCA. We note that attenuation rate of independent scattering is larger than QCA. 
The results show that the extinction rate first increases rapidly with frequency as 
in Rayleigh scattering. Then it departs from Rayleigh scattering around 20 GHz 
and starts to gradually level off. This indicates the importance of including Mie 
scattering effects beyond 20 GHz for (a) = 0.035 cm. Figure 6 compares the 
results of QCA between the extinction rates of monodisperse particles with radii - 
0.035 cm and the case of Rayleigh size distribution with (a) = 0.035 cm. We note 
that the case of Rayleigh size distribution has a considerably larger scattering 
rate than monodisperse case. This indicates the importance of including size 
distribution in modelling scattering from realistic geophysical terrain when grain 
sizes are small. It also means that in laboratory controlled experiments when 
particle sizes are made to be close to monodisperse, the measured scattering can 
be substantially less. The numerical results indicate that the extinction relate 
is at 1.056 cm -1 based on Rayleigh size distribution at frequency =,95 GHz 
meaning that there is is strong multiple scattering of incoherent waves at such 
high frequency even for moderate snow layer thicknesses. Figure 7 shows the 
results of HH backscattering coefficient for a layer thickness of d = 10 cm. It 
shows that distorted Bom backscattering can be as large as 0 dB at 60 GHz with 
(a) = 0 035 cm. The curves shows saturation effects due to the limitations of 
distorted Bom approximation. We expect that if multiple scattering of incoherent 
waves are included, the backscattering should continue to increase above OdB. 
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Figure 5. Coherent wave attenuation rate 2 K{ in cm” 1 as a function 
of frequency. The permittivity of the particles is e s — (3.2 + 
tO.OOl )cq. Rayleigh size distribution is considered with (a) = 
0.035 cm and fractional volume / = 0.3. The cases of indepen- 
dent scattering and QCA are compared. 

QCA (tingle list *=0.035 cm. n Rayleigh lit* <*>=0.035 cm.) 



Figure 6. Coherent wave attenuation rate 2Ki in cm -1 under QCA as a 
function of frequency. The permittivity of the particles is € s = 
(3.2 + i0.001)cq. The fractional volume is / = 0.3. The cases 
of identical sizes of a = 0.035 cm and Rayleigh size distribution 
with (a) — 0.035 cm are compared. 
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Figure 7. Backscattering coefficient &hh in dB of QCA as a function of fre- 
quency. The permittivity of the particles is c$ — (3.2+i0.001)eo- 
The layer thickness is d = 10 cm and incident angle is 30°. The 
cases of identical sizes of a = 0.035 cm and Rayleigh size distri- 
bution with (a) = 0.035 cm are compared. 

In Figs. 8, 9 and 10, the results for the cases of a larger mean grain radius 
of 0.075 cm are shown. The corresponding mode radius is 0.0598 cm. Figure 8 
shows that the extinction rate of independent scattering is still larger than that of 
QCA. However the two curves are much closer than that of Figure 5 showing that 
scattering by larger particles are “more independent.” At frequency = 95 GHz, 
QCA and independent scattering are almost the same. Figure 9 shows that the 
Rayleigh size distribution case is considerably larger than the monodisperse case. 
Strong multiple scattering is to be expected as the extinction rates of QCA for 
Rayleigh size case are large. Figure 10 show the HH backscattering coefficients of 
a snow layer of thickness d — 1 cm. The small layer thickness is chosen to satisfy 
the validity condition of small optical thickness for distorted Born approxima- 
tion. The saturation effect in Fig. 10 is a demonstration of the limitation of the 
distorted Born approximation. Nevertheless, the backscattering coefficients are 
quite large in view of the small layer thickness of d = 1 cm. In Fig. 11, the W 
and HH backscattering coefficients are plotted as a function of incident angles at 
frequency = 37 GHz and d - 10 cm. There is a gradual decrease with incident 
angle. There is little difference between W and HH with the HH return slightly 
larger than W. 
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Figure 8. Coherent wave attenuation rate 2K{ in cm as a function 
of frequency. The permittivity of the particles is e s = (3.2 + 
tO.OOl)e 0 . Rayleigh size distribution is considered with (a) = 
0.075 cm and fractional volume / = 0.3. The cases of indepen- 
dent scattering and QCA are compared. 

qca (single size a*0.075 aad Rayleigh slse <a> * 0.075 cm.) 



Figure 9 . Coherent wave attenuation rate 2 K{ in cm * of QCA as a 
function of frequency. The permittivity of the particles is €$ 
(3.2 + iO.OOl)e 0 . The fractional volume is / = 0.3. The cases 
of identical sizes of a = 0.075 cm and Rayleigh size distribution 
with (a) = 0.075 cm are compared. 
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Figure 10. Backseat tering coefficient in dB of QCA as a function of fre- 

quency. The permittivity of the particles is e 8 = (3.2-M’0.001)6q. 
The layer thickness is d = 1 cm and incident angle is 30°. The 
cases of identical sizes of a — 0.075 cm and Rayleigh size distri- 
bution with (a) = 0.075 cm are compared. 

W and HH retnrn , d=10 cu n Rayleigh six* <a>*0.075 cm. 



Figure 11. Backscattering coefficient a vv and in dB of QCA as a 
function of incident angle. The permittivity of the particles 
is £ $ = (3.2 4* i0. 001)60. The layer thickness is d = 10 cm. 
Frequency is at 37 GHz. The cases of Rayleigh size distribution 
with (a) = 0.075 cm is considered. 
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IV. COMPARISON WITH EXTINCTION MEASUREMENTS OF DRY SNOW BE- 
TWEEN 18 GHZ AND 90 GHZ 

In this section, we apply the dense media QCA-PY-Mie theory to compare with 
extinction measurements of dry snow. Dry snow is a mixture of ice grains and 
air Hallikanen et al. [24] measured the extinction behavior of different types of 
dry snow at four frequencies: 18 GHz, 35 GHz, 60 GHz and 90 GHz. Thus, it is 
important to match the data for a sample over the entire frequency range with 
one set of physical parameters. In this section, we shall match the data of one 
sample (sample number 9 [24]) of large grains. Other samples are being studied 
and comparisons will be reported in the future. 

Recently, grain size distributions of dry snow have been measured extensively 
[27]. Dry snow can be divided into new snow and old snow. New snow has smaller 
mean size and smaller standard deviation while old snow has larger grain size and 
larger standard deviation. The standard deviation of grain size radius can be 
several times larger than the mode radius showing that snow has a broad grain 
size distribution. The maximum grain diameter is between 3 mm to 4 mm. It is 
also observed that the number ratio of large grains to grains with mean size is 
about 1 to 300 [27]. 

In the following, we shall compare with the snow data for sample 9 [24] with 
ground truth of mean grain diameter = 1 mm and snow density = 0.385 gm/cc. 
This corresponds to a fractional volume of 0.423. The reported measurements of 
extinction are 10.2 dB/m, 58.8 dB/m, 247.2 dB/m and 304.0 dB/m respectively 
at frequencies of 18 GHz, 35 GHz, 60 GHz and 90 GHz. Comparisons are made 
with QCA-PY-Mie theory based on two size distributions: the modified gamma 
size distribution and the histogram (Fig. 12) based on information as discussed 
above [27] 

ice grains histogram 

in snow 



0 12 3 


grain diameter (mm) 

Figure 12. The ice grain size histogram in snow in number per c.c. as a 
function of grain diameter. The histogram is used for QCA- 
PY-Mie theory in Fig. 14. There are 15 sizes, with mean grain 
radius of 0.0479 cm and fractional volume of 0.4134. 
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The modified gamma size distribution is 
n(a) = ivja^exp 


with 



(39) 


(40) 


K i = 


3/Q 




(* 2 ) 




(41) 


Thus, the size distributions have four parameters P, Q, K\, and K 2 . How- 
ever, A '2 and A'j can be expressed in terms of mean radius (a) and fractional 
volume / as indicated in (40)-{4l). Thus the alternative four parameters are 
(a),/,P and Q. 

The second size distribution that we used is in Fig. 12 which shows the his- 
togram number density (number per c.c.) for 15 sizes with mean grain radius 
of 0.0479 cm and fractional volume of 0.4134. Note that there are large grains 
between 3 mm to 4 mm with the number ratio to grains of mean size of about 1 


to 300 [27]. 

In Fig. 13, we make a comparison with the data based on the modified gamma 
size distribution with P = 1,Q = 1.4, (a) = 0.05 cm. and / =0.423. The permit- 
tivity of ice grains chosen are €$ = 3.2eo 4- with e^/e o = 0.007, 0.009, 0.011 
and 0.014 respectively at 18 GHz, 35 GHz, 60 GHz and 90 GHz which corresponds 
to a salinity of 0.12 parts per thousand [28]. The results of independent scattering 
are also shown for comparison. The results in figure 13 show that the independent 
scattering results are larger than the data while the QCA-PY-Mie theory shows 
a reasonable match based on the modified gamma size distribution. 

In Fig. 14, we make a comparison with data based on the histogram size distri- 
bution of Fig. 12 and using the same permittivity of ice grains as in Fig. 13. The 
agreement is good. It seems that with size distribution, the QCA-P\ -Mie the- 
ory can explain the frequency dependence of the extinction of dry snow between 
18 GHz and 90 GHz. The results of independent scattering are also shown for 
comparison. Future comparisons will be made when the size distributions of ice 
grains in snow become available [27]. Note that we have used the Percus-Yeyick 
pair functions for non-interacting particles without mutual force. It seems that 
ice grains may have sticking force that may affect the pair distribution functions. 
Besides QCA, the quasicrystalline approximation with coherent potential [6,8] 
should be used. Monte-Carlo simulations of extinction rate can also be made 
based on the solution of Maxwell’s equations [29]. 
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Figure 13. Comparison between QCA-PY-Mie theory and independent 
scattering with snow extinction data at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz of Hallikanen et al. [24] with ground truth 
of mean grain diameter of 1 mm and snow density=0. 385 gm/cc. 
The theory is based on the modified gamma size distribution 
with P = 1, Q = 1.4, < a >= 0.05 cm, and / = 0.423. 
The permittivity of ice grains chosen are e = 3.2«o + with 
e'Z/e o = 0.007,0.009,0.011 and 0.014 respectively at 18 GHz, 
35 GHz, 60 GHz, and 90 GHz which corresponds to a salinity of 
0.12 parts per thousand [28]. 

ground truth of snow: I mm mean 
grain size; 0.385 gm/cc 



Figure 14. Comparison between QCA-PY-Mie theory and independent 
scattering with snow extinction data at 18 GHz, 35 GHz, 
60 GHz, and 90 GHz of Hallikanen et al. [24] with ground truth 
of mean grain diameter of 1 mm and snow density =0.385 gm/cc. 
The theory is based on the histogram of Fig. 12 with < a >= 
0.0479 cm, and / = 0.4134. The permittivity of ice grains cho- 
sen are = 3.2e0 + t£$ with £$/ c 0 = 0.007,0.009,0.011, and 
0.014, respectively, at 18 GHz, 35 GHz, 60 GHz, and 90 GHz 
which corresponds to a salinity of 0.12 parts per thousand [28]. 
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V. CONCLUSIONS 

In this paper, we have illustrated the importance of dense medium theory of 
correlated Mie scattering for studying snow cover at microwave and millimeter 
wave frequencies above 15 GHz. We have also demonstrated the importance of 
including size distributions as there are large difference between the polydisperse 
and the monodisperse cases. The importance of polydispersions has long been 
noted [17], particularly at low frequency when scattering is proportional to particle 
size to the sixth power. Based on size distributions, the dense medium theory is 
able to account for the frequency dependence of extinction of dry snow between 
18 GHz and 90 GHz. The calculations also show that extinction rates are large 
for dense medium so that multiple scattering effects become important even for 
small layer thicknesses. The calculation of the incoherent wave in this paper has 
been limited to that of distorted Born approximation. The subject of multiple 
scattering for the incoherent wave including the effects of correlated Mie scattering 
should be an important future subject for snow cover scattering beyond 15 GHz. 
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